Simple model for the vibrations of embedded elastically cubic nanocrystals 
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The purpose of this work is to calculate the vibrational modes of an elastically anisotropic sphere 
embedded in an isotropic matrix. This has important application to understanding the spectra of 
low-frequency Raman scattering from nanoparticles embedded in a glass matrix. First some low 
frequency vibrational modes of a free cubically elastic sphere are found to be nearly independent of 
one combination of elastic constants. This is then exploited to obtain an isotropic approximation 
for these modes which enables to take into account the surrounding isotropic matrix. This method 
is then used to quantatively explain recent spectra of gold and copper nanocrystals in glasses. 



PACS numbers: 

I. INTRODUCTION 

Low-frequency inelastic light scattering by metallic 
nanoparticles, which is due to their mechanical vibra- 
tions, has been the focus of attention of many researchers 
during the last thirty years^— This scattering is simi- 
lar to surface enhanced Raman scattering from molecules 
close to such a metallic nanostructurei making it an in- 
teresting complementary way to study this complex phe- 
nomenon having many applications for very sensitive de- 
tection. The enhancement due to using laser excitations 
resonant with the dipolar plasmon in such nanostruc- 
tures and the high quality samples available today are re- 
sponsible for low-frequency Raman spectra having an un- 
matched number of features^ compared to non-metallic 
nanoparticles. 

The interpretation of such spectra is very challenging 
as many parameters have to be taken into account at 
the same time. The mechanical vibrations depend on 
the shape of the nanoparticles but also on their inner 
structure and on the surrounding medium. In this work, 
we report on a new approach enabling taking into account 
all these parameters for a spherical nanoparticle having 
a cubic lattice embedded in a glass matrix. 

The vibrational modes of an elastically isotropic sphere 
which is free^ or embedded in an infinite isotropic 
matri:»£i£ are known exactly. When the elastic con- 
stants are not isotropic, a numerical approach such as 
the one known as Resonant UltraSound (RUS)^ can 
find mode frequencies and displacement fields of free 
nanoparticleSfi^— What has been missing up until 
now is a description of the vibrations of an elastically 
anisotropic sphere embedded in an isotropic matrix. The 
case of anisotropic elasticity has been discussed for the 
problem of the scattering of acoustic wave a 13 ' 14 but none 
of these approaches provides the eigendisplacements re- 
quired for modeling the coupling of the vibrations with 
electrons which is at work in all the optical techniques 



used to detect such vibrations. 

Current experimental results such as those presented in 
Ref. have already shown the need for a model without 
any of these limitations. The interpretation in that paper 
took into account the clastic anisotropy to qualitatively 
explain the splitting of the lowest frequency Raman peaks 
but it failed to provide a quantitative description due to 
the significant coupling with the embedding matrix. The 
present work fills that hole by showing that it is possible 
to choose an isotropic approximation of the system for 
the most intense Raman active vibrations which enables 
the prediction of the position of the Raman peaks and 
provides the vibrational displacement fields required for 
the calculation of the Raman intensities. 



II. METHOD 

Exact solutions for the vibrations of isotropic free 
spheres can be classified as spheroidal and torsional and 
will be noted Sf m and T™ m respectively in the following 
with £ and m being the usual angular momentum and 
its z-componcnt and n being an index used to label the 
cigenmodes by increasing frequency starting from n = 1 
as in a previous workJ^ For spheres whose diameter is 
small compared to the wavelength of light, the Raman- 
active vibrations are So and 52 (for every m and 
In the present work, RUS calculations have been used 
to model the vibrations of elastically anisotropic spheres 
as in a previous worki^ by expanding the displacements 
onto x l y 3 z k functions with i + j + k < 20. 

The approach of this paper was inspired by the calcu- 
lations for elastically anisotropic cuboctaedra in Ref. [l6|. 
The idea was to calculate mode frequencies using only 
the speed of sound along a single propagation direction. 
The 5-fold degenerate S\ modes of an isotropic sphere 
are split by cubic elasticity into two degenerate modes 
with Eg symmetry and three with Ti g symmetry (Oh 
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point group) . Due to the symmetry of the displacements 
of the Eg and T 2g vibrations, their frequencies were ap- 
proximated using sound speeds in particular directions 
instead of 3D-averaged ones. It should be noted that 
all the Eg and T 2g vibrations (as well as the A\ g vi- 
brations) are Raman active but only the ones sharing 
a strong similarity with the S 2 modes are expected to 
contribute significantly to the Raman spectra due to the 
surface deformation scattering mechanism.- 

In order to confirm the validity of this simple approach 
for spherical nanoparticles, we use a method similar to 
the one used in a previous work.— Studying the frequency 
changes resulting from a continuous variation of the elas- 
ticity of the material the sphere is made of provides some 
insight into the nature of the vibrations. A cubic ma- 
terial has three independent elastic constants Cn, C\ 2 
and C44 instead of two for isotropic elasticity for which 
C44 = Cll ~ Cl2 , We consider the case of gold nanopar- 
ticles because of available inelastic light scattering ex- 
perimental data to compare with and also because gold 
has a very strong elastic anisotropy making it a good 
system to test the validity of the isotropic approxima- 
tions. We use the following parameters: Cn = 191 GPa, 
C12 = 162 GPa and C44 = 42.4 GPa and mass den- 
sity p = 19.283 g.cm -3 . Since we are mainly interested 
in Raman active vibrations, we will focus on the vibra- 
tions coming from the isotropic S 2 mode which are the 
main features in the low-frequency Raman spectra of gold 
nanoparticles. 

The potential for a great simplification in handling 
these modes, at least in some cases, can be clearly seen 
as follows. We plot eigenfrcqucncies for a gold sphere as 
a function of C44 in Fig. [TJ instead of leaving C44 fixed at 
its normal value for gold. The two lowest frequency E g 
branches, i.e. the fundamental mode and first overtone, 
are very nearly flat. The frequency change for the lowest 
Eg branch is 2% while C44 is multiplied by 4. This in- 
dicates that the corresponding modes approximately do 
not depend on C44 or the associated transverse sound 
speed. This approximate flatness is essential to our ap- 
proach. In practice, it only holds very well for the lowest 
frequency modes, and successively less well for higher fre- 
quency. The projection of the displacements^ of the two 
lowest frequency E g vibrations obtained for yj C44/ p = 
1600 m/s onto those obtained for yf C44/ p = 800 m/s 
is very close to 1 (~ 0.9968). This demonstrates that 
the displacement field of these modes do not significantly 
change with C44 either. Since the mode approximately 
does not depend on C44 , we are free to arbitrarily change 
C44 to a different value which is convenient for us. Specif- 
ically, we can always choose C44 to change the gold into 
an isotropic material, i.e. C44 = 12 . 

The lowest frequency T 2g modes depend on C44 as can 
be seen by their frequency variations in Fig. Q] We also 
consider the variation of their frequencies as a function 
of C\\ — C12 in Fig. [2] while varying either Cn or C\ 2 one 
at a time. There is nearly perfect agreement in this fig- 
ure between varying Cn or C\ 2 . It is also very good for 




FIG. 1: (Color online) Variation of the eigenfrequencies of the 
lowest frequency E g (circles, red online) and Tig (triangles, 
blue online) vibrations of a gold sphere of radius 5 nm when 
varying C44. The vertical arrow indicates the real abscissa for 
gold and the vertical line corresponds to an isotropic system. 



vibrations having other irreducible representations (not 
shown). This demonstrates that Cn — C12 is a good 
choice for a parameter rather than Cn only or C12 only. 
Furthermore, the lowest frequency T 2g vibrations which 
come from the S 2 modes do not depend on the corre- 
sponding transverse sound speed \j Cll 2 ^ ~ ^ as can be 
seen by the almost flat variation. The small variation ob- 
served near the value of gold is due to the anti-crossing 
between the branches of the two lowest frequency T 2g 
branches. While the lowest branch is associated with the 
S 2 mode, the next upper one comes from T3 1 modes. Had 
a less anisotropic material been chosen, this anti-crossing 
pattern would have been less pronounced. Still, as will 
be discussed later, neglecting the mixings between these 
branches is a reasonable choice in many cases. 

In both cases, the mode-dependent isotropic approx- 
imation is obtained by satisfying C44 = Cl1 ~ Cl2 . 
The displacements associated with the mode-dependent 
isotropic approximations can be constructed using sym- 
metry argument a 11 ! 18 from the S 2 , m displacements. S 2 fi 
and S2,2 ~^- 2, ~ 2 are two orthonormal E g modes and 

s ' 2 - 2 J^ 2 '~ 2 , S 2 .i and S 2 ,-i arc the three T 2g orthonor- 
mal modes. As has been shownii, the key parameter for 
the calculation of the S 2 modes is the transverse isotropic 
sound speed although the longitudinal one has also a very 
small contribution. Choosing the isotropic longitudinal 
sound speed is not critical also due to the fact that the 
quasi-longitudinal sound speed does not vary much with 
the propagation direction in gold. As a result, we simply 
used the 3D-averaged longitudinal sound speed as in pre- 
vious works. This value is also the most appropriate one 
for an isotropic approximation of the breathing modes So 
as has been already shown^ 



3 



I 

C5 



250 



200 



150 



100 



50 



t7 



800 900 1000 1100 1200 1300 1400 1500 1600 

((C :: -C 12 )/(2p)) V2 (m/s) 

FIG. 2: (Color online) Variation of the eigenfrequencies of the 
Eg (circles for varying C12 and lines going through the circles 
for varying Cn, red online) and Ti a (triangles for varying 
C12 and lines going through the triangles for varying Cn, 
blue online) vibrations of a gold sphere of radius 5 nm as 
a function of Cn — C12. The vertical arrow indicates the 
real abscissa for gold and the vertical line corresponds to an 
isotropic system. 
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FIG. 3: (Color online) Left: variation of the frequencies us- 
ing RUS (symbols) and Lamb with mode-dependent sound 
speeds (lines) when linearly varying the dj from isotropic 
to anisotropic gold. Radius R = 5 nm. - Right: Matrix 
induced broadening and frequency shift of the same modes. 
vl — a * 5020 m/s, vt — a * 3010 m/s and p = a * 2.97 g/cm 3 
for the matrix and the abscissa a is varied from to 1 so as 
to reach the same value as in Ref. [lq. 



III. APPLICATION 



A. Free gold nanocrystals 



Now the calculation of the eigenfrequencies of the 
modes coming from the spheroidal quadrupolar vibra- 
tions (S2) for a sphere made of a material with cubic elas- 
ticity is separated into two isotropic problems which can 
each be solved exactly. The result of such calculations 
for a gold sphere are presented in the left part of Fig. [3] 
together with RUS calculations for varying anisotropy by 
using Cij(x) = C|?° + x * (Cf™ - C|/°) with < x < 1, 
x being the abscissa, C lso being the isotropic gold stiff- 
ness tensor obtained from the 3D-averagcd sound speeds 
and C am being the anisotropic one. The agreement be- 
tween both kinds of calculations for the lowest three ap- 
proximated branches is very good. The small deviation 
from the lowest frequency T% g branch close to the "free 
anistropic" limit can be attributed to the fact that the 
mixing with the next T 2g branch coming from T3 1 is taken 
into account only in the RUS calculation. This mixing 
has not yet been observed experimentally (it should man- 
ifest as a splitting and intensity sharing between both 
branches). But even in that case, the frequency pro- 
vided by the modified Lamb approach can be seen as a 
good approximation of the position of the expected Ra- 
man peak. As a result, this approach provides a simpler 
description of the vibrations which is quite suitable to 
interpret all the currently available experimental results 
for free nanocrystals. 



B. Matrix embedded nanocrystals 

Since the coupling with an embedding matrix can be 
taken into account for a system having spherical sym- 
metry the isotropic approximations presented before for 
free nanocrystals can be extended to calculate the broad- 
ening and frequency shifts due to a surrounding ma- 
trix. In the following, we focus on calculations using 
the pseudo-modes obtained with the complex frequency 
model (CFM)ii^ Other models such as the core-shell 
model (CSM)i2 would be more suitable for the calcu- 
lation of Raman spectra but won't be considered here as 
we focus only on the positions of the Raman bands. De- 
tails about the validity of the CFM approach using the 
isotropic approximations are presented in Appendix |XJ 

CFM calculations are presented in Fig.[3]and compared 
to the experimental spectra of gold nanocrystals in Fig. 0] 
and Ref. [LQ- Fig. |4] shows the Raman spectra of matrix- 
embedded gold nanocrystals. Details about the sample 
preparation and the spectra acquisition are the same as 
those presented in Ref. [TH with a 64 hour annealing at 
T = 455°C. Fig. @] focuses on the intense lowest frequency 
peak which has been decomposed into a low frequency 
Eg lorentzian (position 141.6 ± 0.7 GHz, full width at 
half maximum (FWHM) 47.5 ± 1.9 GHz) and a higher 
frequency T 2g one (position 184.9 ± 0.4 GHz, FWHM 
63.7 ± 1.2 GHz). 

Fig. [3] (right) presents the evolution from free to 
matrix-embedded gold nanocrystals by varying the pa- 
rameters describing the matrix (mass density and longi- 
tudinal and transverse sound speeds). As discussed be- 
fore, the calculations are expected to be accurate only 



4 



for the lowest E g and T 2g branches but higher frequency 
branches are shown as well for completeness. We used the 
parameters for the matrix which were measured by Bril- 
louin scattering in previous works^^ On the right-hand 
side of Fig. [3l the frequency of the E g branch reaches 
108.3 GHz (FWHM 21.7 GHz) while the T 2g one reaches 
146.2 GHz (FWHM 30.4 GHz). While going from free to 
matrix-embedded nanoparticles, the ratio of the E g and 
T 2g frequencies, which does not depend on the size of 
the nanocrystals, changes from 0.59 to 0.74. Both values 
are in very good agreement with the experimental results 
presented here (ratio 0.77) and in Ref. [H for embedded 
gold nanocrystals^ as well as for free gold nanocrystals 
(ratio 0.62 in Ref. BjJ. This good agreement strongly 
supports the validity of the calculations based on the 
CFM and using the isotropic approximations obtained 
for the free spheres. Moreover, it enables the accurate 
determination of the size of the nanocrystals using the 
inverse proportionality of the eigenfrequencies with the 
diameter of the sphere. 

A similar procedure was used to check the validity of 
our approach with matrix embedded copper nanoparti- 
cles since copper cristallizes in a cubic structure too. Red 
Cu glasses were produced in a similar fashion to the Au 
based glasses described in Ref. HH By annealing near 
the glass transition temperature an initially transparent 
sodo-silicate glass containing minute quantities of both 
CU2O and SnO, nanometric clusters of metallic copper 
are formed. Based on a Maxwell Garnctt description, 
the volume content of Cu under the form of nanoparti- 
cles is estimated as 10 -5 . Similar values are expected for 
the previous gold samples. Low-frequency Raman spec- 
tra were recorded using a laser excitation close to the 
surface plasmon resonance maxima (A = 561 nm for Cu 
instead of A = 532 nm for Au). Due to the observed 
splitting of the lowest frequency band, it is deduced that 
a substantial fraction of the formed Cu nanoparticles are 
monodomain nanocrystals. The density of the embed- 
ding glass is p — 2.43 g.cm~ 3 and the longitudinal and 
tranverse sound velocities are respectively vl = 5800 m/s 
and vt = 3510 m/s, as determined from Brillouin spec- 
troscopy. Note that these matrix parameters differ from 
those of the Au sample (p = 2.97 g.cm -3 , vl = 5020 m/s 
and vt = 3010 m/s). 

The mass density and C^'s for copper were obtained 
from Ref. [20L For free monodomain copper nanocrys- 
tals, the calculated ratio of the E g to T 2g frequencies 
is 0.56 and it increases to 0.66 when embedded in the 
corresponding matrix. The ratio deduced from the fit 
presented in Fig. [5] is 0.71. 

It should be noted that, as clearly demonstrated 
elsewhere^, multiply-twinned particles also contribute to 
the Raman spectra in the same frequency range through 
a broader peak. The broadening of the E g and T 2g 
peaks for matrix embedded nanocrystals prevents the 
clear identification of this additional contribution and 
therefore it is not possible to reliably fit it. However, 
by not taking it into account, the positions of the E g and 
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FIG. 4: (Color online) Low-frequency Raman spectrum for 
matrix embedded gold nanoparticles (crosses). The fit of the 
intense lowest frequency band by two lorentzians and a con- 
stant background is shown with lines. 
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FIG. 5: (Color online) Low-frequency Raman spectrum for 
matrix embedded copper nanoparticles (crosses). The fit of 
the intense lowest frequency band by two lorentzians and a 
constant background is shown with lines. 



T 2g peaks are not very accurate which may explain the 
disagreements between the calculated and fitted ratios. 
For the same reason, the fitted intensities and widths of 
both peaks are seriously affected by the presence of such 
a third contribution. It is therefore important to restrict 
the usage of the fitting procedure used in this work to 
obtain only relatively accurate positions for the E g to 
T 2g peaks. These should of course not depend strongly 
on the exact shape of the peaks (Lorcntzian or Gaussian 
for example). 

Finally, an additional small peak at higher frequency 
can be seen in the spectra presented in Figs U] and [5] 
Within the isotropic approximation, this peak was as- 
signed to the S 2 degenerated modest In the view of 
Fig. El i.e. accounting for the elastic anisotropy and the 
embedding matrix, this peak may tentatively be assigned 
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to either the E g modes deriving from the 5f vibrations 
or the T 2g modes deriving from the 5f modes. Since the 
latter modes are known to have a very small surface de- 
formation, their expected Raman scattering cross-section 
is expected to be negligible. Therefore we assign this 
third small peak to the E g modes deriving from the 5f 
vibrations. 



IV. CONCLUSION 

We have demonstrated that full anisotropic calcula- 
tions are often not required to interpret low-frequency 
Raman spectra from spherical nanoparticles with cubic 
elasticity Instead, simpler isotropic calculations with 
properly chosen sound speeds can provide approximate 
but accurate frequencies for the most intense Raman ac- 
tive vibrations and allow to take into account the medium 
surrounding the nanocrystals. Calculations have been 
successfully compared to experimental results for the case 
of monodomain gold and copper nanoparticles in a glass 
matrix having isotropic elasticity. While the effect of the 
clastic anisotropy on the frequencies of the free vibrations 
was already known, the present work demonstrates the 
measurable impact it also has on the frequency splittings 
and the broadenings of the pseudo-mode for matrix em- 
bedded nanoparticles. Such is important for a reliable 
size evaluation of nanoparticles from their low frequency 
Raman spectra. 
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Appendix A: Validity of the CFM used with the 
isotropic approximations 

Let's consider a gold sphere of radius R . In or- 
der to apply the CFM approach to the case of spher- 



ical nanocrystals having anisotropic elasticity, we have 
to check that the displacement fields originally used in 
the isotropic CFM case are still approximate solutions of 
the wave equation for the anisotropic system for r < Ro 
provided the correct Cy's are used. 

Without any loss of generality, let's first consider the 
Eg modes similar to £2,0 (m = 0). We already know that 
the isotropic approximation presented before is valid for 
the lowest two E g modes. We write the frequencies of 
these modes as v\ = and v 2 = jf- with A < B 
in the following. This isotropic approximation is valid 
whatever the radius of the gold sphere and in particu- 
lar for radii R > Rq. For a given frequency v < 

let's consider two particular spheres of radii Ra = — 
and Rb = 77 so that the first E g mode of the first 
particle and the second E g mode of the second parti- 
cle have the same frequency v. In the isotropic approxi- 
mation, these 5*2,0 displacements are linear combinations 
of two terms: Ul = Vj2(kLr)P2(cos8)e 1,u ' t and ut = 
V x V (r j 2 {kLr)P2 (cos 9)) e zut where jt and Pi are the 
spherical Besscl functions of the first kind and the Leg- 
endre polynomials respectively and u = HlVl = krvr- 
Since two independent linear combinations of the same 
fields (ul and ut) are approximate solutions of the wave 
equation, each field is an approximate solution too for 
r < Rq. Therefore we have checked that it is possible 
to apply the CFM approach in the frequency domain 

v < V\. 

The matrix often increases the pseudo-mode frequen- 
cies up to roughly 50% (i.e. roughly Vl j^ 2 ) in the case 
of a very hard and dense matrix. We want to apply the 
CFM in that frequency range too. The fact that the 
calculations are also valid for R < Rq is of no use here 
because it only demonstrates that the 5*2,0 displacement 
field is a good approximation for the core of the sphere 
but nothing is known close to r = Rq. While we can't ex- 
tend this proof to v > v\ , we can't prove cither that this 
approximation fails very quickly when increasing v above 
v\. For T2 g modes, the situation is also more problematic 
because the branch coming from 5| mixes significantly 
with other T2 g branches coming from 5|, T§ and Sq. 
Despite these limitations, applying the CFM can be use- 
ful at least as a rough estimation of the pseudo-modes 
frequencies. 
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